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Abstract 

We introduce gauge-invariant quark and gluon angular momentum distribu- 
tions after making a generalization of the angular momentum density op- 
erators. From the quark angular momentum distribution, we define the 
gauge-invariant and leading-twist quark orbital angular momentum distribu- 
tion Lg(x). The latter can be extracted from data on the polarized and 
unpolarized quark distributions and the off-forward distribution E(x) in the 
forward limit. We comment upon the evolution equations obeyed by this as 
well as other orbital distributions considered in the literature. 
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It has long been suspected that the quark orbital angular momentum (OAM) plays an 
important role in the structure of the nucleon, even though empirical evidence suggests that 
the constituent quarks are predominantly in the s-wave state in the naive quark model. For 
instance, Sehgal [|IJ noticed, following the work of Ellis and Jaffe [0, that the quark spin falls 
short in accounting for the spin of the proton. He considered seriously the possibility of a 
sizable quark OAM filling in this gap. Ratcliffe studied the angular momentum conservation 
in the parton splitting processes, pointing out that the OAM can be generated by parton 
helicity during the scale evolution ||. Since the EMC data on polarized deep-inelastic 
scattering [[|, a large number of theoretical papers have appeared on the spin structure of 
the proton, some of which have directly and indirectly considered the issue of the quark 
OAM §. 

In Ref. ||, the evolution equation for the quark and gluon OAM was derived in the 
light-cone gauge and coordinates, and in leading-order QCD perturbation theory. Based on 
this, it was concluded that the fraction of the nucleon spin carried by quarks is about h/A 
in the asymptotically large scale limit. Taking into account the polarized DIS data [|],[7|], 
this indicates that approximately h/8 of the nucleon spin resides in the quark OAM in the 
asymptotic limit. In Ref. ||, a proposal was made to gain access experimentally to the 
quark OAM at a finite scale /z 2 by studying virtual Compton scattering in the deep-inelastic 
limit. 

Recently, Hagler and Schafer |§, and Harindranath and Kundu JlO| studied the quark 



OAM distribution as a function of Feynman x in the nucleon. They considered a tower of 
operators, 

^ + (xHd 2 - x 2 id 1 )iD + ■ ■ ■ iD + ip , (1) 

whose matrix elements are used to define the moments of the quark OAM distribution. 
They derived the one-loop evolution equation which has a homogeneous term the same as 



the evolution kernel for the unpolarized quark distribution ||3_ 1|| a generalizing a result found 
earlier for the first moment . The mixing of the quark OAM distribution with the quark 
helicity distribution can also be related to the polarized and unpolarized Altarelli-Parisi 
kernels, as elucidated recently by Teryaev [|T2]]. 

In this paper, we introduce an alternative quark OAM distribution. Our distribution is 
motivated from the quark angular momentum distribution defined from the form factors or 
matrix elements of spin-independent twist-two operators. The new quark OAM distribution 
can be extracted from the polarized and unpolarized quark distributions and the off-forward 
distributions E(x) introduced in Ref. ||. One important feature of the new distribution is 
that it is leading twist and its evolution is completely determined by that of polarized and 
unpolarized quark and gluon distributions to all orders in perturbation theory. 

We start with the familiar spin-independent, twist-two, flavor-singlet (summing over 
quark flavors) quark operators, 

QPm-Hn^ = ^ T 09 _ iD n n ) ^ ^ ( 2 ) 

where all indices /3, fit, n n are symmetric and all traces have been subtracted, as indicated 
by the bracket (■■•)• Note that the n = 1 operator is just the traceless part of the energy- 
momentum tensor of quarks. From O^ 1 ' " At ™, we define a new tower of twist-two operators 
which generalize the angular momentum density, 
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jtf = ^O^ 1 " ^^) - fO^ 1 -^^) - (traces) . 



(3) 



where the subtraction of traces ensures that Jlfwi'-^n belongs to an irreducible represen- 
tation of the Lorentz group. The matrix element of the above operator in the nucleon state 
| -PS), where P and S are respectively the momentum and polarization vectors (P 2 = M 2 , 
S 2 = -M 2 , S-P = 0), is 



PS 



PS 



2J, 



2S p P a 



(^OtfipO pHl . . . pVn _|_ gtoW P CT _p |8] . . . pMn 



9n (n + l)M 2 

+ . . . + e [w^i . . . p/3] _ (traces)) (2vr)V(0) , (4) 



where [a/3] indicates antisymmetrization of the two indices. Here terms with derivatives on 
the 8 function have been omitted because they depend on the center-of-mass motion of the 
nucleon. The matrix element J qn generalizes the nucleon spin carried by quarks J q i. In fact, 
we can define from it a quark angular momentum distribution J q (x) such that 



J q (x)x 



n-1 



dx = J, 



qn • 



(5) 



J q {x) at negative x represents the antiquark contribution. If we limit the support to < 
x < 1, the angular momentum distribution is J q (x) + J q (—x). Obviously one can extend 
all of the above discussions to gluons and define the gluon angular momentum distribution 



Jg(x) 



The matrix element J qn , or equivalently, the angular momentum distribution J q (x) can 
be obtained from the off-forward distribution H(x,£,t) and E(x,£,t) in the forward limit 
£ = t — 0. According to Ref. ||, these distributions are defined in terms of the off-forward 
matrix element of the light-cone string operator, 



i\x 



2tt 



■0 



A 



-71 



" A//2 da n-A(an) , ( A 

w ( - n 



P 



■irr fiu n A 

H q (x,t,t) U(P') flU(P)+E g (x,U) U(P') ^ V U(P) . 



(6) 



where U(P') and U(P) are the Dirac spinors, A = P' - P, P = (P' + P)/2, and t = A 2 . 
Vector n along the light-cone (n 2 = 0) is conjugate to P in the sense that P ■ n = 1. Taking 
the 77,-th moments of the above expression, we get, 



■n, n (P'\0^-^\P) 
H qn (Z,tW(P') iU{P)+E qn {^t)U{P') 



- fTyn ia ^ nfl A v 

2M y ' 



(7) 



where 



H q n^,t) = j^dxx n - l H q (x,^t)dx 
Eqn(£,t) = J dxx n ~ 1 E q (x, £, t)dx 



(8) 
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are the moments of the off-forward distributions. To find their relation with J qn , we write 
down all possible elastic form factors of the twist-two operators after using the constraints 
from Lorentz symmetry and parity and time reversal invariance 

_ [VI _ 

(P'|0^-^|P) = U{P')^ X U{P) J2 A qn , 2i (t)A^ .••A^+ 1 P At2l+2 ...p Mn) 

i=0 

+ U{P') a ^ a U{P) ± B qnj2i (t)^ 2 ■ • • A^P M ■ ..P^ 

i=0 

+ C qn (t) Mod(n + 1, 2) ^U(P')U(P) ■ ■ ■ A' tn) , (9) 

For n > 1, even or odd, there are n + 1 form factors. C n (i) is present only when n is even. 
Substituting the above expression into Eq. (f|), we have 

<V = ^ (/W,o(0) + P g n+1, (0)) . (10) 

Contracting both sides of Eq. (|9]) with n Ml • • ■n fln and comparing the result with Eq. (0), 
we obtain the desired relation 

Jqn= \ (iW (0)0) +^ ? n+l(0,0)) . (11) 

Clearly H q (x, 0, 0) is just the singlet quark distribution q(x). For convenience we abbreviate 
E q (x, 0, 0) as E q {x). Then the above relation can be translated to one for the quark angular 
momentum distribution, 



~x(q(x) + E q (xj) . (12) 
The distribution q(x) has been extracted from high-energy scattering with good accuracy 



J q {x) ^ 



T3| . However, E{x) is unknown at present except for its first moment (see below). 

Since forming spatial moments of an operator as in Eq. (H) does not change its short 
distance behavior, the evolution equation for the angular momentum distributions J q (x) and 
J g (x) is exactly the same as that for the unpolarized quark and gluon distributions, 



d ( Jnq{v) \ = l 7m(" + 1) lqg{ n + !) \ ( J rM 

dln/i 2 I JrM ) V 7 9 ?( n + !) 7 33 (^ + 1) J V Jnfl(/i) 



(13) 



where the anomalous dimension matrix is a perturbation series in a s and the leading order 
result can be found in Ref. ||11|| . Of course, according to Eq. (|12"|) the same conclusion follows 



from the evolution of the off-forward parton distributions ||. 

Now we can define the quark orbital angular momentum (OAM) distribution by looking 
at the structure of M q y/3fll '" fln . Since it is already an irreducible operator under Lorentz 
transformations, the generalized OAM operator that it contains is not irreducible and hence 
the nucleon matrix element of the latter cannot be a Lorentz scalar. This is an intrinsic 
limitation to the significance of the notion of the orbital angular momentum in a relativistic 
quantum theory. However, if we are limited to a class of coordinates in which the nucleon 
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has a definite helicity, we still can define the orbital angular momentum and its distribution 
in Feynman x [Q. In the following discussion, we assume the nucleon is moving in the z di- 
rection and is polarized in the helicity +1/2 state. We consider only the specific components 
a = 1, (3 = 2, \i\ = • • • = fi n = + of the M q operator and take its spatial integral, 



/ ^Mf+-+ = / d^ (f 1 ^) - e 2 M( 



i+-+) 



(14) 



where we have introduced the light-cone coordinates: ^ = (£°±£ 3 ) / y/2. The spatial integral 
element is <i 3 £ = d 2 £±d£~ . One can perform partial integrations to make all derivatives act 
on the right, 
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n + 1 J 1 

+ ip>y + iD 2 ■ ••iD + ip + 



■ ■ + $^ + iD~ 



iD 2 ^ -(1^2). (15) 



The first term in the bracket can be further manipulated, 

1 



V>7 iD + ■ ■ ■ iD + ip 



n 



%l)(^ 2 iD + - 1 + iD 2 )iD + ■ ■ -iD+tp + 



+ ipiD + ■ ■■iD + { 1 2 iD + - 7+iD 2 )^ 



+ 



1 r T 



7? 



+ ^ + iD + ■■■iD 2 i) 



(16) 



The second term in brackets on the right hand side is the same as the remaining terms in 
Eq. (|i~5|). So we concentrate on the first term in the above equation. Using the following 
identity, 



^ 2 iD + - 



z^7 + 7 2 - 7 + 7 2 i J p 



(17) 



and partially integrate the first derivative term, we have, 

d 3 £ \$(y 2 iD + - >y + iD 2 )iD + ■ ■ ■ iD + i) + •■■ + $iD+ ■ ■ ■%D + (y 2 iD + - y + iD 2 )^ 



n 



-%- I £^7 7 7 ■ ■ ■ iD + ^ + / d^ ipy + (x l Y)(igF p+ )y p iD + ■ ■ ■ iD + ^p 



+ 



+ (n - l)tpy + iD + • ■ ■ iD + (x 1 ^ 2 ){igF p+ ) lp i) 



(18) 



where we have used the quark equations of motion to eliminate some terms. 
Putting all the pieces together, we have 



dW 



(19) 



where, 



n 



V3 

d 3 £ip~f + —iD- 



- / d 3 £ Uy + (xHD 2 - x 2 iD 1 )iD~ 
n J L 



iD + )^ 
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AL+-+ 



+ ••• + ^ + iD + ■■■iD + (x 1 iD 2 - x 2 iD 1 )ip 
1 



n(n + 1) 
+ ••• 



d 3 £ Vn^OV - x 2 l l ) (igF p+ -f p )iD + ■ ■ ■ iD + ip 
+ %[) 1 + iD + ■ ■ ■ iD + {x 1 ^ 2 - arV) {igF p+ lp )^ . 



(20) 



The matrix element of S q h is clearly related to the polarized quark distribution Aq(x) 
which has been the major focus of the polarized deep-inelastic scattering experiments in the 



last ten years 



seems to be a natural gauge-invariant generalization of the OAM 



operator. However, from its definition, it is clear that L+"' + contains not only the leading 
twist but also the higher-twist contributions. In particular, under a change of renormaliza- 
tion scale, it mixes with AL+"' + . Only the sum + AL+"' + evolves as a leading-twist 

operator. Therefore, we define the leading- twist, gauge-invariant generalization of the OAM 
operator as 



T+-- 



Its matrix element in the nucleon state is 

2L 



(PS\L+- + \PS) 



qn 



71 + 1 



AL + " H 



2P+...P+(2tt)V(0) 



(21) 



(22) 



Let us emphasize again that the simple structure of the matrix element is only possible in 
the specific class of frames that we defined earlier. From the above, the OAM distribution 
function L q (x) can be defined as, 



Lq{x^x dx 



(23) 



which can be expressed in terms of the quark angular momentum and helicity distributions, 



L q (x) = J q (x) - -Aq(x) . 



(24) 



The definition of the matrix elements of h and its relation to the quark helicity distri- 
bution is standard, 



(PS\St- + \PS) 



n + 1 



2P+...P+(2tt) 3 5 3 (0) 



J dxx^Aqix) = S n . 



Combining Eqs. (|T2,|24), we have 



L q {x) 



x 



(q{x) + E q 



x 



Aq(x) 



(25) 



(26) 



Thus, the quark OAM distribution can be determined entirely from the polarized and un- 
polarized singlet quark distribution and the off-forward distribution E(x). 

The evolution equation for L qn (fi) is straightforward. The evolution of S n is well-known, 
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dlia.fi 2 n 



A^ qq (n)H n + Aj qg (n)Ag Ti 



(27) 



where Ag n are the moments of the gluon helicity distribution and the anomalous dimensions 
can be found in Ref. II] . Combining the above with Eq. (|T3"D, we have the evolution equation 
for L n 



■'qni 



+%g( n + l ) J gn ~ ^Aj qg (n)Ag n , 



(28) 



which is again valid to all orders in perturbation theory. 

In view of the above discussion, we now have a new perspective on the evolution equations 
obtained by Hagler and Schafer, and Harindranath and Kundu. The quark OAM operators 
defined in these works (we will call it Ljj~ " + ) differ from our "' + by terms which depend 
on transverse polarizations of the gluon potentials (A 1 , A 2 ). Our calculation indicates that 
these extra terms do not mix in the L q0 " + operator at the one-loop level, although at 
higher-loops they must do so that L+" 



maintain an all-order simple evolution in Eq. (f28|). 
Therefore at the one-loop level, L q0 " + and L^"' + have an identical homogeneous part in 
their evolution equation. Taking into account the extra factor of in our definition of 
the OAM matrix elements, we deduce that L° qn + AT, n /(n + 1) has the same homogeneous 
evolution as J, 
that f0 

these, we immediately deduce that 



qn at one-loop level. Similar analysis for the gluon angular momentum shows 
Lg n + 2Ag n / {n + 1) has the same homogeneous evolution as J gn at one loop. From 



d 



din /j, 2 \ L u g M 
a s (jj) 1 



tts(/f) 

2vr 



7°> + l 



, o 



+ 



n 



2n n + 1 



lqq\ 
"fgq( n + 1) 




iqq 
2A 7g Vr/ 



AE n (//) 

A^n(^) 



+ 



(29) 



at one loop. The explicit terms agree with the known result RIC]. The ellipses denote the 
matrix elements of operators with explicit transverse gluon fields which are present even at 
the one-loop level. Because of these operators, the evolution equation is not guaranteed to 
have a similar simple structure at higher loops. 

To determine the quark OAM distribution, one has to measure E(x) from off-forward 
hard scattering processes. Deeply-virtual Compton scattering has been proposed as a way 
to measure it [H. High-energy diffractive vector- meson production can also provide access 
to the distribution [0. The first moment of E(x) is related to the anomalous magnetic 
moment, 



J E{x)dx = k u + K d + K s 



(30) 



The experimental data on this is 0.33 ±0.39, a fairly small number 16|| . This could indicate 



a small size of E(x). If E(x) is indeed small, L q {x) can be determined entirely in terms of 
Aq(x) and q(x). However, there is no sound theoretical reason to neglect E(x) at present. 
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